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Abstract 

We define a manifold M where objects c 6 M are curves, which we parame- 
terize as c : S 1 — > R" (71 > 2, S 1 is the circle). Given a curve c, we define the 
tangent space T C M of M at c including in it all deformations h : S 1 — > R n of c. 

In this paper we study geometries on the manifold of curves, provided by 
Sobolev-type metrics H J . We study H 3 type metrics for the cases j = 1, 2; we 
prove estimates, and characterize the completion of the space of smooth curves. 

As a bonus, we prove that the Frechet distance of curves (see [MM06b]) co- 
incides with the distance induced by the "Finsler L°° metric" defined in §2.2 in 
[YM04b]. 



1 Introduction 

Suppose that c is an immersed curve c : S 1 — ► R™, where S 1 C Ft 2 is the circle; we 
want to define a geometry on M, the space of all such immersions c. 

The tangent space T C M of M at c contains all the deformations h £ T C M of the 
curve c, that are all the vector fields along c. Then, an infinitesimal deformation of the 
curve c in "direction" h will yield (on first order) the curve c(u) + eh(u). 

For the sake of simplicity, we postpone details of the definitions (in particular on 
the regularity of c and h, and the topology on M) to §2. 

We would like to define a Riemannian metric on the manifold M of immersed 
curves: this means that, given two deformations h,k G T C M, we want to define a 
scalar product (h, k) c , possibly dependent on c. The Riemannian metric would then 
entail a distance <1{cq,c\) between the curves in M, defined as the infimum of the 
length Len(7) of all smooth paths 7 : [0, 1] — ► M connecting cq to c\. We call 
minimal geodesic a path providing the minimum of Len(7) in the class of 7 with fixed 
endpoints. (1) 

At the same time, we would like to consider the curves as "geometric objects"; to 
this end, we will define the space of geometrical curves B=M/Diff(S' 1 ), that is the 
space of immersed curves up to reparametrization. To this end, we will ask that the 
metric (and all the energies) defined on M be independent of the parameterization of 
the curves. 

B and AI are the Shape Spaces that are studied in this paper. 



A number of methods have been proposed in Shape Analysis to define distances 
between shapes, averages of shapes and optimal morphings between shapes. At the 
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same time, there has been much previous work in Shape Optimization, for example Im- 
age Segmentation via Active Contours, 3D Stereo Reconstruction via Deformable Sur- 
faces; in these later methods, many authors have defined Energy Functionals E(c) on 
curves (or on surfaces), whose minima represent the desired Segmentation/Reconstruction; 
and then utilized the Calculus of Variations to derive curve evolutions to search minima 
of E{c), often referring to these evolutions as Gradient Flows. The reference to these 
flows as gradient flows implies a certain Riemannian metric on the space of curves; but 
this fact has been largely overlooked. We call this metric H°, and define it by 

(h,k) H0 = T I (h(s),k(s))ds 



L 



s 1 



where h, k g T C M, L is the length of c, ds=|c(#)| d6 is integration by arc-parameter, 
and (h(s), k(s)) is the usual Euclidian scalar product in IR™ (that sometimes we will 
also write as h(s) ■ k(s)). 

For example, the well known Geometric Heat Flow is often referred as the gradient 
flow for length: we show how it is indeed the gradient flow w.r.t. the H° metric. 



Example 1 Let c be an immersed curve, and h be a deformation of c. Let the differen- 



tial operator D s =-irdg be "the derivative with respect to arclength". Let 



len(c)*/ |c(0)| 
Js 1 



de (i) 



be the length of the curve; we recall that 
dlen(c) 



dh 



[ (D s h ■ T) ds = - I (h- Die) ds 
Js 1 Js 1 



(2) 



where T — D s c is the tangent to the curve, D^c is the curvature, intended as a vector. 

Let C = C{6,t) be an evolving family of curves trying to minimize L: we re- 
call moreover that the resulting geometric heat flow ( also known as motion by mean 
curvature) 



dt 

is well defined only for positive times. 

By comparing the above flow to the definition ofH°, we realize that this flow is the 
gradient descent (up to a conformal term 1 / len(c) ): 

dC 1 „ , , . 

a^-i^) v " olcn(c) 

If one wishes to have a consistent view of the geometry of the space of curves in 
both Shape Optimization and Shape Analysis, then one should use the H metric when 
computing distances, averages and morphs between shapes. 

Surprisingly, H° does not yield a well define metric structure, since the associated 
distance is identically zero (2) . 

Moreover, some simple Shape Optimization tasks are ill-defined when using the 
H° metric: 



* 2 'This striking fact was first described in [Mum]; it is generalized to spaces of submanifolds in [MM05] 
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Remark 2 Let g : R™ — > R fc , to 



))=t/ = £/ 



(Tzere again L=\en(c)); then 



daVg[9{c)) = jf V 9 (c)h + 9 (c)(D s h.T)ds- 
L Js 1 

^ j j(c) ds j (D B h-T) ds = 



dh 

L 



= -/ Vg(c)h+(g(c)-avg(g(c)))(D s h-T)ds (3) 

s 1 



If the curve is in the plane, that is n = 2, then we define the normal vector N _L T 
by rotating T counterclockwise, and define scalar curvature k so that D^c = kN; 
then, integrating by parts, the above becomes 

davg{g{c)) 1 f dg ,dg >. 

-(g(c)-avg(g(c)))(h-D 2 s c) ds = 
= \ J (|f (c)N - «( 5 (c) - avg 07(c)) )) (h- N)d 



s 



Suppose now that we have a Shape Optimization functional E including a term of 
the form avg (g(c)); let C = C{6, t) be an evolving family of curves trying to minimize 
E; this flow would contain a term of the form 

^ = ...(g(c{s))-g) K N... 

unfortunately the above flow is ill defined: it is a negative-time heat flow on roughly 
half of the curve. We present two simple examples. 

• If for example g(x) = x, then avg (g(c)) = avg (c) is the center of mass of the 
curve. Let us fix a target point v £ R 2 . Let E(c)=^\avg (c) — v\ 2 be afunctional 
that penalizes the distance from the center of mass to v; by direct computation 

dE . . . davg (c) 

— = { avg{c)-v).—^ r = 

1 



L 



((c-v) ■ (iV — k(c — c))){h ■ N) ds 



(where c = avg (c) for simplicity); so we conclude that the H° gradient descent 
flow is 

dC 

— = -V H oE(c) = ({v - c) ■ N)N - kN({c -c)-(v- c)) (4) 



(up to a part tangent to the curve). 
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Let P={w : ((«; — c) ■ (v — c)) > 0} be the halfplane 
that is the region of the plane that is "on the v side " w.r.t. 
c. This gradient descent flow (4) does move the center 
of mass towards the point v: indeed there is a first term 
((v — c) • N)N that moves the whole curve towards v; 
and a second term that tries to decrease the curve length 
out of P and increase the curve length in P: and this is 
ill posed. 

• Similarly if g{x) = \x — c\ 2 , then E(c)=avg (|c — c| 2 ) is the standard deviation 
of the curve. The derivative is 

\ f (<(c -c)-N)- n(g(c) - avg (g(c)) )) (h ■ TV) ds 

The flow to minimize this should be 

dC 

— = -<(c - c) ■ N)N + K N(g(c) - avg (g(c)) ) 

and this is ill posed where the curve is inside of the circle of center avg (c) and 
radius y^avg (5(c)) 

The above phenomenon is also visible in many applications, where the Active Con- 
tour curve would "fractalize" in an attempt to minimize the task energy. For this reason, 
a regularization term is often added to the energy: this remedy, though, does change 
the energy, and ends up solving a different problem. 




§l.i Sobolev-type Riemannian Metrics 

To overcome this limitation, in [SYM05a] and [SYM06] we proposed a family of So- 
bolev-type Riemannian Metrics 

Definition 3 Let c £ M, L be the length of c, and h,k € T C M. Let A > 0. We assume 
h and k are parameterized by the arclength parameter of c. We define, for j > 1 
integer, 

i) . (h, k) H3 * (h, k) H0 + XL 2j (Dih, D{k) H0 

ii) . (h, k) S) d d aV g (h) ■ avg (k) + \L 2] (Dih, D{k) R0 

where again avg (h) =j- J sl h(s) ds and D 3 S is the j-th derivative with respect to ar- 
clength. (3) 

It is easy to verify that the above definitions are inner products. Note that we have 
introduced length dependent scale factors so that the above inner products (and corre- 
sponding norms) are independent of curve rescaling. 

Changing the metric will change the gradient and thus the gradient descent flow; 
this change will alter the topology in the space of curves, but the change of topology 
does not affect the energy to be minimized, or its global minima; whereas it may reg- 
ularize the flows, and avoid that the flows be trapped in local minima; many examples 
and applications can be found in the survey [SYM05b]. 

* 3) Note that (h, k) H o = avg (h ■ k) so the difference in the two metrics is in using avg (h ■ k) instead of 
avg (h) ■ avg (k) 



§l.ii Previous work 
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§l.ii Previous work 

In [YM04a, YM04b, YM05] we addressed the problem of denning a metric in the 
space M of parameterized immersed curves c : S 1 — ► R" (with special attention to 
the case n = 2); we discuss the general setting of Finsler and Riemannian metrics, 
and related problematics; we propose a set of goal properties, we discuss some models 
available in the literature. Eventually we proposed and study conformal metrics such 
as 

(h,k) H o = len(c) J (h(s),k(s))ds (5) 

we prove results regarding this metric, and in particular, that the associate distance is 
non degenerate. We also proved that, supposing that only unit length curves with an 
upper bound on curvature are allowed, then there exist minimal geodesies. 

The same approach was proposed independently by J. Shah in [Sha05], who more- 
over proved that in the simplest case given by (5), the only minimal geodesies are given 
by a curve evolution whose velocity is proportional to the curve normal vector field. 

The problem has been addressed by Michor and Mumford in [MM06b], who pro- 
pose the metric 

(h,k) H o A =[ (1 + AK 2 (s))(h(s),k(s))ds (6) 

where k is the curvature of c, and A > is a fixed constant; they prove many results 
regarding this metric; in particular, that the associate distance is non degenerate, and 
that completion of smooth curves is contained between the space Lip of rectifiable 
curves, and the space B V 2 of rectifiable curves whose curvature is a bounded measure. 

More recently in [SYM05a] and [SYM06] we studied the family of Sobolev-type 
metrics defined in 3. 

In [SYM05a] we have experimentally shown that Sobolev flows are smooth in the 
space of curves, are not as dependent on local image information as H° flows, are 
global motions which deform locally after moving globally, and do not require deriva- 
tives of the curve to be defined for region-based and edge-based energies. In general, 
Sobolev gradients can be expressed in terms of the traditional H° gradient, that is, we 
have the formulas 

V H ^E = V H oE*K X!n (7) 
\7 Rn E = V H oE*k x , n (8) 

for suitable convolutional kernels K\^ n , K\ iTl . We have moreover shown mathemati- 
cally that the Sobolev-type gradients regularize the flows of well known energies, by 
reducing the degree of the P.D.E. 

Example 4 For example, in the case of the elastic energy E(c) = J k 2 ds = J \D 2 c\ 2 ds, 

the H° gradient is \7 H oE = LD S (2D^' c + 3\D 2 c\ 2 D s c) that includes fourth order 
derivatives; whereas the H 1 -gradient is 

Lc ss + L{\D ss c\ 2 D s c) * K XA (9) 

that is an integro-differential second order P.D.E. 

In [SYM06] we have shown numerical experiments on real-life cases, and shown 
that the regularizing properties may be explained in the Fourier domain: indeed, if we 
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calculate Sobolev gradients Vh^ E of an arbitrary energy E in the frequency domain, 
then 

^«> = t^Sf f " iEZ ,10) 

and 

^~E(0) = ^E(O), V^E(l) = forZeZ\{0}, (11) 

(see eqn. 13 for the precise definition of Fourier coefficients). It is clear from the 
previous expressions that high frequency components of Vh°E{6) are increasingly less 
pronounced in the various forms of the H n gradients. 

A family of metrics similar to the first example above (but for the length dependent 
scale factors) is currently studied in [MM06a]: the Sobolev-type weak Riemannian 
metric on Imm(iS , R 2 ) 

(h,k) Gn = [ y^(D%Dik)ds= [ (A n {h),k)ds where 

n 

A n (h) = A n , c (h) = J2(-V i Ds i (h) 

i=0 

in that paper the geodesic equation, horizontality, conserved momenta, lower and upper 
bounds on the induced distance, and scalar curvatures are computed. Note that this 
norm is equivalent to the norms in 3: see remark 19. 



Many other approaches to Shape Analysis are present in the literature; for example, 
much earlier than the above, Younes in [You98] had proposed a computable definition 
of distance of curves, modeled on elastic curves. 

A different approach to the study of shapes is obtained when the shape is defined 
to be a curve up to reparametrization, rotation, translation, scaling. 

For example Mio, Srivastava et al in [KSMJ03, MS04] use a different choice of 
curve representation: they represent a planar curve c by a pair of angle-velocity func- 
tions c(u) = cxp(0(u) + iO(u)) (identifying R 2 = (C), and then defining a metric 
on (<fi, 8). They propose models of spaces of curves where the metrics involve higher 
order derivatives in [KSMJ03]. See the proof of thm. 24 for an example comparison of 
the two approaches. 



2 Spaces of curves 

As anticipated in the introduction, we want to define a geometry on M, the space of all 
immersions c : S 1 — > R™. 

We will sometimes distinguish exactly what M is, choosing between the space 
Imm(S' 1 , R") of immersions, the space Imm/(S' 1 , R") of free immersions, and Emb(S' 1 , R") 
of embeddings (see §2.4,§2.5 in [MM06b]). 

We will equip M with a topology r stronger than the C 1 topology: then any such 
choice M is an open subset of the vector space C 1 (S 1 , R") (that is a Banach space), 
so it is a manifold. 

The tangent space T C M of M at c contains vector fields h : S 1 — > R" along c. 
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Note that we represent both curves c G M and deformations h G T C M as functions 
S 1 — ► R n ; this is a special structure that is not usually present in abstract manifolds: 
so we can easily define "charts" for M: 

Remark 5 (Charts in AI) Given a curve c, there is a neighbourhood U c ofO G T C M 
such that for h G U c , the curve c + h is still immersed; then this map h i— > c + h is the 
simplest natural candidate to be a chart of $ c : U c — > M; indeed, if we pick another 
curve c G AI and the corresponding ?7g such that J7g H U c =/= 0, then the equality 
<5> c (h) = c + h = c + h = ^^(h) can be solved for h to obtain h = (c — c) + h. 

The above is trivial but is worth remarking for two reasons: it stresses that the topology 
r must be strong enough to mantain immersions; and is a basis block to what we will 
do in the space B^t defined below. 

We look mainly for metrics in the space M that are independent on the parameter- 
ization of the curves c: to this end, we define these spaces of geometrical curves 

Bi = B^S^TR 2 ) = Imm(S ,1 ,R 2 )/Diff(S' 1 ) 

and 

Bij = Bij{S\lR 2 ) = Imm / (S' 1 ,R 2 )/Diff(S' 1 ) 

that are the quotients of the spaces Imm/ and Imm/(S' 1 , R 2 ) by Dif^S 1 ); alternatively 
we may quotient by Diff + (S' 1 ) (the space of orientation preserving automorphisms of 
S 1 ), and obtain spaces of geometrical oriented curves. 

Remark 6 (on model spaces and properties) We have two possible choices in mind 
for the topology r to put on M: the Frechet space of C°° functions; or a Hilbert space 
such as standard Sobolev space iP (S 1 — > R"). 

Suppose we define on M a Riemannian metric: we would like Bi to have a nice 
geometrical structure; we would like our Riemannian Geometry to satisfy some useful 
properties. 

Unfortunately, this currently seems an antinomy. 

If M is modeled on a Hilbert space H 3 , then most of the usual calculus carries 
on; for example, the exponenential map would be locally a diffeomorphism; but the 
quotient space M/Diff[S ) is not a smooth bundle, (since the tangent to the orbit 
contains c and this is in H°~ l in general!). 

If AI is modeled on the Frechet space of C°° functions, then the quotient space 
M/Dijf(S ) is a smooth bundle; but some of the usual calculus fails: the Cauchy 
theorem does not hold in general; and the exponenential is not locally surjective. 

Suppose in the following that r is the Frechet space of C°° functions, for simplicity; 
then Bi j is a manifold, the base of a principal fiber bundle while Bi is not (see in §2.4.3 
in [MM06b]). 

To define charts on this manifold, we imitate what was done for M: 

Proposition 7 (Charts in B^/) Let II be the projection from Immf(S 1 ,lR) to the 
quotient Bi j. 

Let [c] G Bi j: we pick a curve c such that 11(c) = [c]. We represent the tangent 
space T^Bij as the space of all k : S 1 — > R™ such that k(s) is orthogonal to c(s). 

Again we can define a simple natural chart <E>j c ] by projecting the chart <E> C (defined 
in 5): the chart is 

$ [c] (£0=n( c (.) + fc(-)) 
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that is, it moves c(u) in direction k{u); and it is easily seen that the chart does not 
depend on the choice of c such that 11(c) = [c]. We can solve $r c i (fc) = <j>^] (fc) ( this is 
not so easy to prove: see [MM06b], or 4.4.7 and 4.6.6 in [Ham82]). 

Define a Finsler metric F on M; this is a lower semi continuous function such that 
F(c, •) is a norm on T C M. 

If 7 : [0,1] — > M is a path connecting two curves cq, c\, then we may define a 
homotopy C : S 1 x [0, 1] — > R™ associated to 7 by (7(0, = j(v)(6), and viceversa. 

Definition 8 (standard distance) Given a metric F in M, we could consequently de- 
fine the standard distance of two curves cq, c\ as the infimum of the length 

f F( 7 (i),7(*))di 
Jo 

in the class of all 7 connecting Co, C\. 

This is not, though, the most interesting distance for applications: we are indeed inter- 
ested in studying metrics and distances in the quotient space B=M/Diff(S' 1 ). 
We suppose that 

Definition 9 the metric F(c, h) is "curve-wise parameterization invariant", that is, it 
does not depend on the parameterization of the curves c 

then F may be projected to B=M /~D\S(S X ); we will say that F is a geometrical 
metric. 

Consider two geometrical curves [cq], [ci] E B, and a path 
7 : [0, 1] — > Bi connecting [cq], [c\\: then we may lift it to a 
• ; ; • cH>' homotopy C : S 1 x [0, 1] — * in this case, the homotopy 
■' ; ; ■' will connect a cq o </>o to ci o X , with 0o, 0i G Diff(S' 1 ). 

Since i 71 does not depend on the parameterization, we can 
; factor out <fra from the definition of the projected length. 

To summarize, we define the 

Definition 10 (geometric distance) Given Co, ci, we define the class A of homotopies 
C connecting the curve cq to a reparamete rization c\ o <j> of the curve C\, that is, 
C(u,0) = cq(u) and C(u,l) = Ci(<p(u)). We define the geometric distance dp of 
[ c o] 5 [ c i] m B'^AI/Diff(S 1 ) as the infimum of the length 

Len F (C)= f F (C(-, v), d v C(-,v)) dv 
Jo 

in the class of all such C £ A. ^ 

Any homotopy that achives the minimum o/Lerii?(C) is called a geodesic. 

We call such distances dp(co, Cx), dropping the square brackets for simplicity. (6) 
We provide an interesting example of the above ideas. 

(4) If we use Diff + (S 1 ) to define B then <f> must be orientation preserving as well. 

(5) Note the difference between Len(C) and len(c), that was defined in eqn. (1). 

* 6) We are abusing notation: these dp are not, properly speaking, distances in the space M, since the 
distance between c and a reparameterization c o tj> is zero. 




§2.i L°° -type Finsler metric and Frechet distance 



9 



§2.i L 00 -type Finsler metric and Frechet distance 

We digress from the main theme of the paper to prove a result that will be used in the 
following. 

For any fixed immersed curve c and e 5 1 , we define for convenience 7Tjv : R" — > 
R n to be the projection on the space N(9) orthogonal to the tangent vector D s c(9), 

Tr N{g)W = w-(w,D s c(6))D s c(6) VweR". (12) 

Consider two immersed curves Co and c\\ the Frechet distance df (as found in 
[MM06b])is defined by 

Definition 11 (Frechet distance) 

d/(c ,ci)=inf sup \ci((j>(u)) - c (u)\ 

<t> u 

where u € S 1 and (j) is chosen in the class of diffeomorphisms of S . 

This is a well defined distance in the space Bi (that is not, though, complete w.r.t. this 
distance: its completion is the space of Frechet curves). 

Another similar distance was defined in §2.2 in [YM04b] by a different approach, 
using a Finsler metric: 

Definition 12 (Finsler L°° metric) If we wish to define a norm F(c, •) on T C M that 
is modeled on the norm of the Banach space L°°(S l — ► R n ), we define 

F°°( c , h)Z\\«sh\\L~ = sup W N{e) h(9)\ 
e 

We define the distance doo(co, Ci) as in 10. 

Section §2.2.1 in [YM04b] discusses the relationship between the distance and the 
Hausdorff distance of compact sets; we instead discuss the relationship between df and 
doo : indeed we prove that df — doo . 

Theorem 13 df = d^. 

Proof. Fix c and c\, and define A as in 10. 

We recall that d^ is also equal to the infimum of 

f 1 dC 
rfoo(co,ci)= inf / sup I — (9,v)\ dv 

as well (the proof follows immediatly from prop. 3.10 in [YM04b]) 

Consider a homotopy C = C(u, v) G A connecting the curve Co to a reparameter- 
ization c\ o cf) of the curve c\ : 

sup \ci(cf>(u)) ~ co(u)\ = sup|C(u, 1) - C(u,0)\ = 

U U 

nl dC 



sup 



f dc. « _ r ac. . 
J ^ u ^ dv ^ J s 7 ~dv- {u ' v) 



dv 
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On the other side, let 

C^(0,tO=(l -«)<*>(*) 
be the linear interpolation: then 

dC* . 

(that does not depend on v) so that 

.1 dc< p 



sup 



(it, v) dv 



ci(0(w)) - 


co(u) 




ac* 


= / sup 

J It 


dv M 



dv 



and then, for that particular homotopy C*, 

Lenoo (C*) = sup | ci (4>(u) ) - c (u) | 

u 

we compute the infimum of all possible choices of <fr and get that 

doc(c ,ci) = inf Lcn 00 (C) < infLen 00 (C*) = inf sup |ci(^(u))-c (u)| = d/(co,Ci) 



c 



□ 



The theorem holds as well if use orientation preserving diffeomorphism Diff + (S' 1 ) 
both in the definition of the Frechet distance and in the definition of L°°. 



3 Sobolev-type W metrics 

We start by generalizing the definition in 3. Fix A > 0. Suppose that h G L 2 , then we 
can express it in Fourier series: 



h(s) = 2^ hQ) ex P 



lez 



2m, 

T 



-Is 



(13) 



where h £ i? 2 (Z — > C); and similarly for k. 

For any a > 0, given the Fourier coefficients h, k : Z — > C of h, k, we define the 
fractional Sobolev inner product 



(M>*a=X;(27rQ a «fc(0-*(0 



(14) 



lez 



that is independent of curve scaling; then we can define 

(h, k) Ha =avg {h-k) + A (/i, k) Hg 

(h, k) ^ a =avg (h) ■ avg (fc) + A (h, k) Ha 



(15) 



When a = j integer, these definition coincide with the one in 3. So, for any a > 0, 
we represent the Sobolev-type metrics by 



(h, k) Ha = + X{2nl) 2a )h{l) ■ k(l) 

lez 

(h, k) fj a = h(0) ■ k(0) + H2nl) 2a h(l) ■ k(l). 



(16) 
(17) 



lez 
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Remark 14 Unfortunately for j that is not an integer, the inner products (therefore, 
norms) are not local, that is, they cannot be written as integrals of derivatives of the 
curves. An interesting representation is by kernel convolution: given r 6 R + , we can 
represent them, for j integer j > r + 1 /4, as 



(h,k)ffr= J J D°h{s)K(s - §)D 3 k($)dsds 

that is, {h 1 k) ^ r = (D 3 ' h, K * D^k} „ Q , for a specific kernel K, ( * denotes convolu- 
tion in S 1 w.r.t. arc parameter). 

Remark 15 The norm \\h\\jjj has an interesting interpretation in connection with ap- 
plications in Computer Vision. 

Consider a deformation h £ T C M and write it as h — avg (h) + h: this decomposes 

T C M = R n D C M (18) 

with 



D c M=^h : S 1 -> R" | avg (h) = o} 



If we assign to R™ its usual euclidean norm, and to D C M the scale-invariant Hq 
norm, then we are naturally lead to decompose as in eqn. (15), that is 

\\h\\% a = \avg{h)\^+X\\h\\ 2 HS (19) 

This means that the two spaces IR™ and D C M are orthogonal w.r.t. H a . 

In the above, R™ is akin to be the space of translations and D C M the space of non- 
translating deformations. That labeling is not rigorous, though! since the subspace of 
T C M that does not move the center of mass avg (c) is not D C M, but rather 

\h: j h+ (c-avg{c))(D s h-T)ds = oj 

according to eqn. (3). 

The decomposition (19) is at the base of a two step alternating algorithm for min- 
imization of H 1 Sobolev Active Contours (see $3.4 in [SYM05a]), where the tracking 
of contours is done by, alternatively, minimizing an energy on curves translations, and 
then on curves deformations in D C AI with the metric \\h\\ H i. The resulting algorithm 
is independent of the choice of X. 

Note that {h, h) Ha is a norm on D C M (by (21)), and it is a seminorm and not a 
norm on T C M. 

We define norms as 



F Hi (c,h) = \\h\\ Hj = y/{h,h) Hi , F 6j (c,h) = \\hWfjj = \f(Jhh) 
and consequently we define distances d^i and as explained in 10. 



Hi 
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§3.i Stuff 

We improve a result from [SYM05a] <7) : we show that the norms associated with the 
inner products H ■? and W are equivalent. We first prove 

Lemma 16 (Poincare inequalities) Pick h : [0, L] — > R", weakly differentiable, with 
h(0) = h(L) (so h is periodically extensible) so that 

h{u)-h{{))= ti(s)ds = - ti(s)ds 

JO Ju 

then derive these equations 

h(u)-h(Q) = \{^f n h'(s)ds- j h'(s)d^ => 
=> avg(h)-h(0) = 



\avg(h)-h(Q)\ < ±£ ^£ \h'(s)\ds + £\h'(s)\dsj du 
1 f L ( f L \ . 1 ' L 



2L 



o \Jo 



\h'(s)\ds \ du = - / \h'(s)\ds 



so that (by extending h and replacing with an arbitrary point) 

<-L 



awp\h(u) -avg(h) | < \ [ \ti(s)\ds (20) 
u 2 Jo 



'o 

the constant 1/2 is optimal and is approximated by a family of h such that h'(s) = 
o(l[ 0je) (s) - l[ e> 2e)(s)) when e — * (for a fixed a e IR"). (8) 

By wsmg Holder inequality we can then derive many useful Poincare inequalities 
of the form \\h — avg (K) || p < c Pj1JJ || h'\\ q . By Fourier transform we can also prove for 
p = q = 2 that 

\h( S )-avg(h)\ 2 ds<j^-j o |^)( S )| 2 d S (21) 

where the constant C2,2.j = {L/2ir) 2: ' is optimal and is achieved by h(s) = as\n(2irs/ L) 
(with a S R"j. (9) 



Proposition 17 Fix a smooth immersed curve c : S 1 — ► R n , /ef L = len(c). By 
Holder's inequality, we have that \avg (h) | 2 < J Q L |/i(s)| 2 ds 50 f/iaf < 
On the other hand, 

\ £ \h(s) - avg (h) | 2 ds = i ^ |/i( S )| 2 d S - |avg (ft) | 2 (22) 

(7) and we provide a better version that unfortunately was prepared too late for the printed version of 
[SYM05a] 

®1j4 (x) is the characteristic function, taking value 1 for x £ A,0 for x £ A. 

* 9) This h is not the only solution; for n = 2 we also have h(s) = (cos(27rs/L), sin(27rs/L)). 
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so f/iaf (by the Poincare inequality (2\)) 



\h\\h 



L 

< \avg 
Consequently, 



y\h(s)\ 2 + XL^-^h^is)] 2 ds 

Lj 

\ [ L \h( S )-avg(h)\ 2 ds+ ( L \L 2 3- l \h^(s)\ 2 ds+\avg{h){ 



c)i 2+i2j - i ((^ +i )ri" ra wi 2d ^ i ^?'"'i 



2 

Hi 



/l + (27r) 2 iA 



More in general 



Proposition 18 For i = 0, . . . , j, choose ao > anrf a; > wtt ao + ao > and 
> 0. Define a H 3 -type Riemannian norm * 10 ' 



3 pL 

avg(h)\ 2 + y2a i L M - 1 \h^{s)\ 2 ds (23) 

„-_n JO 



»=0 



f/ien a/Z smc/i norms are equivalent. 

Moreover, choose r with 1 < r < j, and choose bo > 0, 6, > vw'fn 60 + 60 > 
0, 6 r > 0: fnen fne norm || h\\ ( a yj is stronger than the norm \\h\\^ r . 

Proof. The proof is just an application of (22) and of (21) (repeatedly); note also that 
for 1 < i < j equation (21) becomes 

\h^(s)\ 2 ds< (27r)2j _ 2i y \h®(s)\ 2 ds (24) 

since avg =0. □ 

So our definitions of || • \\ HJ and || • || ^ are in a sense the simpler choices of a 
Sobolev type norm that are scale invariant; in particular, 

Remark 19 the H 3 type metric 

\\h\\ 2 M M fj2\h^(s)\ 2 ds 

studied in [MM06a ] is equivalent to our choices, 

bi\\-\\ HJ <||-||m<6 2 ||.|| Hj 
but the constants b±, 6 2 depend on the length of the curve. 

Following this proposition, we will prove some properties of the H 1 metric, and 
we will know that they can be extended to H 1 and to more general iP-type metrics 
defined as in (23). 



We prove this fundamental inequality (25): 

''"'the scalar product can be easily inferred 
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Proposition 20 Fix a smooth immersed curve c : S 1 — > IR", let L = len(c). 
We rewrite for convenience 

\\h\\ 2 H1 > XL 2 (h',h') Htt = XL \h'(s)\ 2 ds = X J \c(u)\du J \ti [u)\ 2 \c{u)\ du 

where hi = D s h; then by Cauchy- Schwartz 

\c(u)\du [ \h'(u)\ 2 \c(u)\du > ( [ \h'(u)\\c(u)\du 



Suppose now that C(u, v) is a smooth homotopy of immersed curves C(-, v): then set 
h(u, v) = d v C(u, v) so that D s h = D s d v C = d uv C /\d u C\. Summarizing the above 

RCM)|| ffl > / \d uv C(u,v)\du (25) 



As argued in 18, the above result extends to all H° -type norms (23) . 

We related the iJ x -type metric to the L°° type metrics 

Proposition 21 The H 1 metric is stronger than the L°° metric defined in 12. 

As a consequence, by theorem 17 and 13, the H 3 and H J distances are lower 
bounded by the Frechet distance (with appropriate constants depending on X). 

Proof. Indeed, by (20) there follows 

sup \n N(fi) h(0)\ < sup \h(8)\ < |avg (h) \ + \ [ \ti\ds < 
ft a ^ ./ 



< |avg (h) | + J \h'\Hs < V2J |avg (h) P + I J \h'\Hs 



(7Tjv was defined in eqn. (12)). For example, choosing A = 1/4, 

Foo(c,h) < V2\\h\\Hi 



□ 



We also establish relationship between the length len(c) of a curve and the Sobolev 
metrics: 

Proposition 22 Suppose again thatC(u, v) is a smooth homothopy of immersed curves, 
let L(y)= len(C(-, v)) be the length at time v; then 

d v L = J(d uv C,-^)du< J\d uv C\du<-L\\C v (;v)\\ m 

by (25). 

We have many interesting consequences: 

|L(l)-L(0)|<-^Len(C) (26) 

where the length Lcn(C) of the homotopy/path C is computed using either H 1 
or H 1 (or using any metric as in (23) above, but in this case the constant in (26) 
would change). 
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• Define the length functional c i— * len(c) on our space of curves; embed the space 
of curves with a H 1 metric; then the length functional is Lipschitz. 

• The "zero curves" are the constant curves (that have zero length); these are points 
in the space of curves where the space of curves is, in a sense, singular; by the 
above, the "zero curves" are a closed set in the H 1 space of curves, and an 
immersed curve c is distant at least len(c)vA from the "zero curves". 

But the most interesting consequence is that 

Theorem 23 (Completion of B\ w.r.t. H 1 ) let d H i be the distance induced by H 1 ; 
the metric completion of the space of curves is contained in the space of all rectifiable 
curves. 

Proof. This statement is a bit fuzzy: indeed is not a distance on M, whereas in 
B objects are not functions, but classes of functions. So it must be intended "up to 
reparametrization of curves", as follows. 

Let (c„)„ e iN be a Cauchy sequence. Since does not depend on parametriza- 
tion, we assume that all c„ are parametrized by arc parameter, that is, |<9ec„| = l n 
constant in 8. By proposition 21 , all curves are contained in a bounded region; 
since lcn(c„) = 2wl n by proposition 22 above, the sequence l n is bounded. So the 
(reparametrized) family (c n ) is equibounded and equilipschitz: by Ascoli-Arzela the- 
orem, up to a subsequence, we obtain that c n converges uniformly to a Lipschitz curve 
c, and |<9ec| < lim„ l n . □ 

We also prove that 

Theorem 24 Any rectifiable planar curve is approximable by smooth curves according 
to the distance induced by H 1 . 

Proof. Let c be a rectifiable curve, and assume that it is non-constant. 

As a first step, we assume that c is not flat, that is, the image of c is not contained 
in a line in the plane. We sketch how we can approximate c by smooth curves. The 
precise arguments are in section 2.1.4 in [YM04b]; see in particular the proofs of 2.12 
and 2.15 in the appendix. Since the metric is independent of rescaling, we rescale c, 
and assume that \dec\ = 1. 

We identify S 1 with [0, 2tt). Let in the following L 2 = L 2 ([0, 2tt\. We define the 
measurable angle function to be a function r : [0, 2tt) — > [0, 2ir) such that dgc(0) = 
(cos t (6), sin t {9)). We define 



(this is similar to what was done in Srivastava et al. works on "Shape Representation 
using Direction Functions", see [KSMJ03]). 

As proved in 2.12 in [YM04b], S is a manifold near r. As shown in the proof 
of 2.15 in [YM04b], there exists a function ir : V — > S defined in a neighbourhood 

'"'The concept is clarified by introducing the concept of horizontality in M, that we must unfortunately 
skip for sake of brevity 



S= {ra 2 ([0,2,]) |0(r) 



(0,0)} 



where : L 2 — > R 2 is defined by 
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V C L 2 of t such that, if f(s) 6 L 2 is smooth, then 7r(/)(s) is smooth. Let /„ be a 
smooth approximation of r, with /„ — > r in L 2 ; then g n =^(fn) — * 7". Let then 

G n (fl,t)=7r(t7- + (l-% n )(fl) 

be a the projection on 5 of the linear path connecting r to g n . Since 5 is smooth in 
V, then the L 2 distance ||r — g n \\ is equivalent to the geodesic induced distance; in 
particular, 

lim E s G n = 



where 



E S G= f \\d t G(;t)\\hdt 
Jo 



is the action of the path G in S C L 2 . 

The above G n can be associated to an homotopy by definining 

C„(M)=c(0)+ / (cos(G„(0,i)),sin(G„(M)))d0 



note that G„(s, 0) = c(s) and C n (s, 1) is a smooth closed curve. 
We now compute the H 1 action of C n , 

Efii(G„)= C \\d t C n \\ 2 Hi &t= [ [ |9 t G„| 2 + |^ s 9 t C„| 2 d s di 
Jo Jo Jo 

Since any G n (-, t) is by arc parameter, then D s dtC n = d s tC n so 

D s d t C n = N(s)d t G n (s,t) 

where 

N(s)=(- sin(G„(s, t)), cos(G„( S , t))) 

is the normal to the curve; so the second term in the action E H i (G n ) is exactly equal 
to Es(G n ), that is, 

pi pi plu 

Effi(G„)= / ||ftC„||^i dt = / / |9 t G„| 2 d s di + ^ s (G„) 
Jo Jo Jo 

We can also prove that J Q J Q 27r |9tG„ | 2 ds dt — > 0, so Ejyi (G n ) — > 0, and then 

limLcn^fi (G„) = 

n 

As a second step, to conclude, we assume that c is flat, that is, the image of c is 
contained in a line in the plane; then, up to translation and rotation, 

c(0) = ( Cl (0),O) 

since c is by arc parameter, t\ = ±1. Let then / : [0, 2-7r] be smooth and with support 
in [1, 3] and /(2) = 1; let moreover 

c(e,tp( Cl (e),tf(e)) 

so 

\d e c\ = y/i + (f(e)) 2 > i 

and then we can easily prove that 

hen H i (G) < oo 

moreover, any curve C(-,t) for t > is not flat, so it can be approximated by smooth 
curves □ 



§3.ii The completion of M according to H 2 distance 
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§3.ii The completion of M according to H 2 distance 

Let d(co, ci) be the geometric distance induced by H 2 on M (as defined in 10). Let 
E{c)= J \D 2 c\ 2 ds be defined on non-constant smooth curves. We prove that 

Theorem 25 E is locally Lipschitz in M w.r.t. d, and the local Lipschitz constant 
depends on the length of c. 

As a corollary, all non-constant curves in the completion of C°°{S 1 — > IR") ac- 
cording to the metric H 2 admit curvature as a measurable function, and the curvature 
satisfies E{c) < oo. 

Viceversa, any non-constant curve admitting curvature in a weak sense and satis- 
fying E(c) < oo is approximable by smooth curves. 

The rest of this section is devoted to proving the above three statements. 
Fix a curve Co; let Lq= len Co be its length. 

By eqn. (26) and prop. 18, we know that the "length function" c i— > len(c) is 
Lipschitz in M w.r.t the distance d, that is, 

| len Co — lenci| < a\d(cQ, ci) 

where a\ is a positive constant (dependent on A). 
Choose any c\ with d(co, ci) < Lo/(4ai). 

Let C(9, t) be a time varying smooth homotopy connecting cq to (a reparametriza- 
tion of) ci; choose it so that Len C < 2d(co, Ci); then Len C < Lq/(2cli). 

Let L(t)= len C(-, t) be the length of the curve at time t. Since at all times t G 
[0, 1], d(co,C(-,t)) < L Q /{2a 1 ), then \L(t) - L \ < oii /(2ai) = L /2; in particu- 
lar, 

L /2 < L(t) < L 3/2 . 

By using this last inequality, we are allowed to discard L(t) in most of the following 
estimates. 



We can ||/||^/ = 



ds and 



N(tp\\D 2 s d t C(-,t)\\ = ^ J \D*d t C\ 2 d6 
for convenience; using this notation, we recall that 



\\d t C\\ H . = ^\L{tfN(t) 2 + —\\d t C\\ 2 ■ 

so \\d t C\\ H 2 > V\L 3 / 2 N{t). 

Up to reparametrization in the t parameter, we can suppose that the path t C(-,t) 
in M is by (approximate) arc parameter, that is ||(9 t C|j#2 is (almost) constant in t; so 
we assume, with no loss of generality, that ||i9 f C||^f2 < 2d(co, c\) for alH £ [0, 1], and 
then N(t) < a 2 d(c ,ci) where a 2 = 2/ v /(i /2) 3 A. 

We want to prove that 

E(ci) - E(c ) < a 5 d(c ,ci) 



where the constant a$ will depend on Lq and A. 
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By direct computation 

d t E(C(-,t))= J \D 2 C\ 2 (D s d t C,D s C) ds + 2 J(D 2 C,d t D 2 C) ds 

we deal with the two addenda in this way: 
i). by Poincare inequality (20) we deduce 



sup\D s d t C\ <±f \D 2 s d t C\ds < J \D- 2 s d t C\ 2 ds = sfW)N{t) 



since avg (D s dtC) = 0. 

So we estimate the first term as 



\D 2 C\ 2 (D s d t C,D s C) ds < E(C)y/L^)N(t) 



ii). The commutator of D s and dt is (D s dtc, D s c)D s : indeed 



1 i « , fa 1 \a n a (d t d e c,d e c) 

— — dgd t + (d t ——)d e = D s d t — - — i 

\d e c\ \d g c\ |<9 e c| J 

D s d t ~(D s d t c, D s c)D s 



so 

d t D 2 C = D s dtD s C - (DsdtC, D S C)D 2 C = 

= D 2 s d t C - D s ((D s d t C, D S C)D S C) - (DsdtC, D S C)D 2 C 
= D 2 s d t C - ((D 2 s d t C, D S C)D S C) - ((DsdtC, D 2 S C)D S C) - 
-2((D s dtC,D s C)D 2 s C) 

so (since |D S C| = 1) 

\\d t D 2 C\\ < 2\\D 2 d t C\\ + 3\\D 2 C\\sup\Dsd t C\ 

that yields an estimate of the second term 



(D 2 C, d t D 2 C) ds < y/E\C)(2N(t) + 2^E(C)^L{t)N(t) 



by using Cauchy-Schwartz. 
Summing up 



\d t E(C(-,t))\<2^E(C)N(t) + 'mC)VW)N(t) 
or, since y/x < 1 + x, 

\d t E(C(-, t))\ < 2N(t) + 2E(C)N(t) + 3E(C)^/I{t)N(t) 
We recall that N(t) < a 2 d(co, c\), L(t) < Lq2>/2, so we rewrite the above as 
\d t E(C(-,t))\ < 2a 2 d(c Q , ci) + 2E{C)a 2 d(c , a) + 3E(C) ai a 2 d(c , a) 
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with 04 = J Lq3/2. Apply Gronwall's Lemma to obtain 

E{ci) < (e(co) + 2ci2d(co, Ci)J cxp ^(2 + 3a i )a 2 d(c , ci)J . 

Let 

9(y)=(E( c o) + 2a22/) cxp H2 + 3a 4 )a 2 y) 

then E{c\) < g(d(co, Ci)); since g is convex, and g(0) — E(cq), then there exists a 
a$ > such that g(y) < E(cq) + a$y when < y < Lo/(4ai); since we assumed 
that d(co, Ci) < Lo/(4ai), then 

E(ci) < E(c ) + a 5 d(c , Ci) . 

Note that 05 is ultimately dependent on Lq and A. 
This ends the proof of the first statement of 25. 

To prove the second statement, let (c n )n£N be a Cauchy sequence. Since djji < 
adj{2, then as in the proof of 23, we assume that, up to reparametrization and a choice 
of subsequence, c n converges uniformly to a Lipschitz curve c. 

Let Lq — len c. We have assumed in the statement that c is non-constant; then 
L > 0. 

Again, the "length function" c ^ len(c) is Lipschitz, so we know that the sequence 
len(c„) is Cauchy in R, so it converges; moreover the "length function" c 1— > len(c) is 
lower semicontinuous w.r.t. uniform convergence, so lim„ len(c„) > len(c) > 0. So 
we assume, up to a subsequence, that 2Lq > len(c„) > Lq 

We proved above that, in a neighbourhood of c of size Lo/(8ai), the function 
E(c)= J \D 2 s c\ 2 ds is Lipschitz; so we know that the sequence E(c n ) is bounded, and 
then (since curves are by arc parameter and len(c„) > Lq) the energy J \dgc\ 2 ds is 
bounded: then dec n are uniformly Holder continuous, so by Ascoli-Arzela compact- 
ness theorems, up to a subsequence, dgc n (9) converges. 

As a corollary we obtain that lim n lcn(c n ) = lcn(c), that c is parametrized by arc 
parameter, and that D s c n (6) converges to D s c(8). 

Since the functional J \dgC n \ 2 ds is bounded in n, then by a theorem in [Bre86], c 
admits weak derivative dgC and J \d@c\ 2 ds < 00, and equivalently, J \D 2 c\ 2 ds < 00. 

For the third statement, viceversa, let c be a rectifiable curve, and assume that it is 
non-constant, and E(c) < 00. Since the metric is independent of rescaling, we rescale 
c, and assume that |9 u c(u)| = 1. 

We express in Fouries series 

c(u) = l n exp(inu) (27) 

nGZ 

(by equating S 1 = R/27r), then we decide that 

C(u, t)= l n exp(mn — f(n)t) (28) 

nGZ 

with f(n) = f(—n) > and lim/(n)/log(n) = 00; (for example, f(n) = \n\ or 
f(n) = (log(|n| + 2)) 2 ): then C(-,t) is smooth for any t > 

We want to prove that, for t small, C(-,t) is near c in the H 2 metric; to this end, let 
C be the linear interpolator 

C(u, t, r)=(l - t)c{u) + tC(u, t) = J2 Le inu (l - r + re"^"'*) (29) 

nez 
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we will prove that 



J (J \d T C\ 2 ds + XL 4 J \D 2 s d T C\ 2 ds^ dr < S(t) (30) 



where lim t ^o S(t) = 0, and L is the length of C(-, t, r). 
We need some preliminary results: 

• we prove that 

/ \d uu c~d uu C\ 2 du<S 1 (t) 
Js 1 

where lim t _»o <M*) = 0, uniformly in r e [0, 1]; we write 



(31) 



and since 



/ \duuC- d uu C\ 2 du = 27TT 2 V |Z„| 2 M 4 (1 - e- /(n)t ) 

nez 

27rV|/„| 2 |n| 4 = S(c)= / |d u „c| 2 ds<oc 
Z7i Js 1 



and lim t _> (l - erf^) 2 = 0, we can apply Lebesgue dominated convergence 
theorem. 

• We prove that 

\d u c-d u G\<8 2 (t) (32) 
where lim^o S^it) = 0, uniformly in u and re [0,1]; indeed 

\d u c - d u C\ <r^2 IU»I(1 - e" /(n) *) < 



y ngZ y n6Z : n^0 

again we apply Lebesgue dominated convergence theorem. 

• By the above we also obtain that for t small, 

3/2 > \d u C\ > 1/2 uniformly in r, u (33) 

• We can similarly prove that 

\c-C\<6 3 (t) (34) 

By direct computation 

j-)2g q 9 uur C (d uu C , d u C)d UT C 

s T ~ \d u c\ 2 

but then, for t small, by (33), 

\D 2 s d T C\ < 4\d UUT C\ + 2A\d uu C\\d UT C\ 
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We use the fact that 

duurC = d uu C — d uu c , d UT C = d u C — d u c , d T C = C — c, 
so by eqn. (31) and eqn. (32) 

\D%C\ 2 dsdT < o x (*i(t) + E{c)8 2 {t)) 



and by eqn. (34) J J \d T C\ 2 dsdr < 8Ja(f). Eventually we combine all above to 
bound eqn. (30) by setting S(t) = a 2 S 3 (t) + Xa 2 {5i(t) + E(c)S 2 (t)). 
This concludes the proof. 
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